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Abstract 

The paper intends to offer a general overview on what the concept of integrabihty means for 
a nonhnear dynamical system and how the symmetry method can be applied for approaching 
I it. After a general part where key problems as direct and indirect symmetry method or optimal 

■ system of solutions are tackled out, in the second part of the lecture two concrete models of 
^ I nonlinear dynamical systems are effectively studied in order to illustrate how the procedure is 

working out. The two models are the 2D Ricci flow model coming from the general relativity 
and the 2D convective-diffusion equation . Part of the results, especially concerning the optimal 
\^ • systems of solutions, are new ones. 

Keywords: Lie symmetries, invariants, similarity reduction. 

^! 1 Integrability and symmetries. Key aspects. 

d . 1.1 The concept of integrability for dynamical systems 

Dynamical systems described by nonlinear partial differential equations are frequently used to model 
a wide variety of phenomena in physics, chemistry, biology and other fields pj. The modelling 
. process includes to find solutions of those partial differential equations. If these solutions exist, the 

r — \ diffential system is said to be integrable. Sometime it is difficult to find a complete set of solutions 

■ and it is quite enough if one can decide on the integrability of the system. There are many methods 
which can be used to fulfill this aim: the Hirota's bilinear method, the Backlund transformation 
method, the inverse scattering method, the Lax pair operator, the Painleve analysis and others 
[2]. Each method has its own significant properties. For example, while the Lax and the Painleve 
methods are moreover testing the integrability, the Hirota's bilinear method is very efficient for 
the effective determination of the multiple soliton solutions for a wide class of nonlinear evolution 
equations [3]. As a conclusion, to decide that a nonlinear differential equation is integrable, one of 

^ ■ the following situation shouls appear: 

[i) the existence of a number of functionally independent first integrals/invariants equal to the 
order of the system in general and half that for a Lagrangian system as a consequence of Liouville's 
Theorem; 

{ii) the existence of a sufficient number of Lie symmetries to reduce the partial differential 
equation to an ordinary differential equation; 

[in) the possession of the Painleve property [1]. 
In this lecture the first two criteria will be investigated. 

1.2 The symmetry method for solving dynamical systems 

Many natural phenomena are described by a system of nonlinear partial differential equations (pdes) 
which is often difficult to be solved analytically, as there is no a general theory for completely solving 
of the nonlinear pdes. One of the most useful techniques for finding exact solutions of the dynamical 
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systems described by nonlinear pdes is the symmetry method. On the one hand, one can consider 
symmetry reduction of differential equations and thus obtain classes of exact solutions. On the 
other hand, by definition, a symmetry transforms solutions into solutions, and thus symmetries 
can be used to generate new solutions from known ones. 

Initially the symmetry method for solving partial differential equations was developed for what 
is currently known as the Lie (classical) symmetry method (CSM). We shall present now a short 
introduction to this approach |5]. 

Let us consider a n-th order partial differential system: 

A,(x,^.(")[x]) = (1) 

where x = {x*, i = \,p} <Z BP represent the independent variables, while u = {u", a = l^q} <Z W 
the dependent ones. The notation ti^") designates the set of variables which includes u and the 
partial derivatives of u up to n-th order. 

The general infinitesimal symmetry operator has the form: 

The n-th extension of ([2]) is given by: 

U^-^ = U + j2Y.ct>i{xM'^))^ (3) 

a=l J J 

where 

am Q 

u') = f4) 

Also, in dH) the second summation refers to all the multi-indices J = (ji, ...jm),with 1 < < 
p,l < m < n .The coefficient functions 0^ are given by the following formula: 



1=1 1=1 



in which 



"^'^ dx^ dx'^dx^^dx^'^ ..dx^" 

d"" 



(7) 



Vj - P,,P,,...D,„ - ^—^-—^ (8) 

The Lie symmetries represent the set of all the infinitesimal transformations which keep invariant 
the differential system. The invariance condition is: 

f/^"^[A] iA=o=0 (9) 

The characteristic equations associated to general symmetry generator ([2]) have the form: 

dx^ dx'P dv} du'^ , , 

C C (Al (Pq 

By integrating the characteristic system of ordinary differential equations (|1U|) . the in- 
variants Ir, r = 1, (p + g — 1) of the analyzed system can be found. They are identified with the 
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constants of integration. Following this way, the set of similarity variables is found in terms of 
which the original evolutionary equation with p independent variables and q dependent ones can 
be reduced to a set of differential equations with {p + q — 1) variables. These are the similarity 
reduced equation which generate the similarity solution of the analyzed model. 

There have been several generalizations of the Lie symmetry method which include: 

1) the non-classical symmetry method (NSM) (also referred to as the conditional method) of 
Bluman and Cole [6], 

2) the direct method of Clarkson and Kruskal [llj, 

3) the differential constraint approach of Olver and Rosenau [12j 

4) the generalized conditional symmetry method due to Fokas, Liu and Zhdanov [13j . 

The basic idea of the nonclassical method is that should be augmented with the invariance 
surface condition: 

g-(x,n«) = 0jx,n)-^f(:E,u)^=O, a=l,q (11) 

1=1 

The g— tuple Q = {Q^ , Q^, ■■■Q'^) is known as the characteristic of the symmetry operator ([2|). The 
invariance condition ([9]) must be applied taking into account that the constraints (jlip do exist. 
The number of determining equations for the infinitesimals (,^{x,u), (/)„(x,n), appearing in the 
nonclassical method is smaller than for the classical method. The main difficulty of this approach 
is that the determining equations are no longer linear. On the other hand, the NSM may produce 
more solutions than the CSM, since any classical symmetry is a nonclassical one, but not conversely. 

The direct method represents a direct, algorithmic, and nongroup theoretic method for finding 
symmetry reductions. The relationship between this direct method and the nonclassical method 
has been discussed in many papers (e.g., 0, [8]). In particular, Levi and Winternitz [9] established, 
using a group-theoretic explanation, that all new solutions obtained by the direct method can also 
be obtained by the nonclassical method. In fact, it has been shown in [1^ that the similarity 
solutions corresponding to the nonclassical groups should in general constitute a larger family than 
that obtained by the direct method. 

The differential constraint approach proposed a generalization of the nonclassical method. Its 
promoters shown that many known reduction methods, including the classical and nonclassical 
methods, partial invariance, and separation of variables can be placed into a general framework. In 
their formulation, the original system of partial differential equations can be enlarged by appending 
additional differential constraints (side conditions), such that the resulting overdetermined system 
of partial differential equations satisfying compatibility conditions. 

As well, in further efforts to find new symmetries of PDEs which would lead to additional 
new invariant solutions, much work has been done in the area of higher-order symmetries. In 
particular, for an evolution equation in two independent variables and one dependent variable has 
been introduced in [13] the method of generalized conditional symmetries ( GCS) or conditional 
Lie-Backlund symmetries. 

1.3 Optimal system of solutions 

In general, when a differential equation admits a Lie group Qr and its Lie algebra Cr is of dimension 
r )1, one desires to minimize the search for invariant solutions by finding the nonequivalent branches 
of solutions. This leads to the concept of optimal system. 

It is well known that for one-dimensional subalgebras, the problem of finding an optimal sys- 
tem of subalgebras is essentially the same as the problem of classifying the orbits of the adjoint 
transformations. 

In Ovsiannikov [T3], the global matrix of the adjoint transformations is used in constructing the 
one-dimensional optimal system. 
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In Olver [5], a slightly different technique is employed: it consists in constructing a table, named 
the adjoint table, showing the separate adjoint actions of each element in Cr as it acts on all the 
other elements . 

The procedure reported in Ruggieri and Valenti [15], is a mixed of the above procedures and 
consists in constructing the global matrix of the adjoint transformations by means of the adjoint 
table. 

One of the advantages of the symmetry analysis is the possibility to find solutions of the original 
pdes by solving odes. These odes, called reduced equations, are obtained by introducing suitable 
new variables, determined as invariant functions with respect to the infinitesimal generators. 

On the basis of the infinitesimal generators of the optimal systems of the Lie algebras of analyzed 
model, we can construct the reduced odes of the given model and find exact solutions. 



1.4 The inverse Lie symmetry problem 

Usually, the direct symmetry problem of evolutionary equations is considered for finding their exact 
solutions. It also known as the classical symmetry method. Firstly, it consists in determining the 
Lie symmetry group corresponding to a given evolutionary equation. Then, using the characteristic 
equations could be obtained the Lie invariants associated to each symmetry operator. Further these 
invariants, following the reduced similarity procedure, determine the reduced equation which could 
be solve and generates the similarity solution of the analyzed model. 

Also, the inverse symmetry problem yL6j could be made. We ask the question: what is the largest 
class of evolutionary equations which are equivalent from the point of view of their symmetries?. 
So, this problem could be solved by imposing a concrete symmetry group to a general analyzed 
model. With this condition, the general symmetry determining equations could be solved and allow 
to determine all concrete models which admit the same Lie symmetry group. 

Let us consider a 2D dynamical system described by a second order partial derivative equation 
of the general form: 

Ut = A{x,y,t,u)uxy + B{x,y,t,u)uxUy + C{x,y,t,u)u2x + D{x,y,t,u)u2y + 

+E{x, y, t, u)uy + F{x, y, t, u)ux + G{x, y, t, u) (12) 

with A{x,y,t,u), B{x,y,t,u), C{x,y,t,u), D{x,y,t,u), E{x,y,t,u), F{x,y,t,u), G{x,y,t,u) arbi- 
trary functions of their arguments. 

The general expression of the Lie symmetry operator which leaves (|12p invariant is: 

d d d d 

U{x,y,t,u) = (p{x,y,t,u)— +i{x,y,t,u)— + r]{x,y,t,u)— + (t){x,y,t,u)— (13) 

Through loss the generality we can choose in the previous expression (p = l.Then, the generator 
(fT3|) becomes: 

d d d d 

U{x,y,t,u) = — + i{x,y,t,u)— + r]{x,y,t,u)— + (p{x,y,t,u)— (14) 

Following the symmetry theory [5], the second extension C/^^^ of ()13p has to be considered and the 
invariance condition of the equation ()12p is given by the relation: 

= C/^^) [ut - A{x, y, t, u)uxy - B{x, y, t, u)uxUy - C{x, y, t, u)u2x - D{x, y, t, u)u2y - 

-E{x, y, t, u)uy - F{x, y, t, u)ux - G{x, y, t, u)] (15) 
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The previous relation has the equivalent expression: 

= -AfUxy - BtUxUy - CtU2x - DtU2y - EfUy - FfU^ - Gt - A^^u^y - B^iu^Uy - 

-Cxiu2x - Dxiu2y - Exiuy - Fxiux - Gxi - Ay-quxy - ByrjUxUy - Cyr]U2x - DyriU2y - 

-EyTJUy - FyTJUx " G yl] - Au(t>Uxy " Bu4>UxUy " Gu(t>U2x " Du4>U2y - Eu(f)Uy - Fu<f>Ux - 

-Gu<l) + Ac^'^y - Get?"" - Dcf^y - Bc^'^uy - FcP"^ - B^yux - E^y (16) 
The functions cf)^ , cf)^ , 4>y , 4>'^^ , (f)'^y , (j)^y will be determined using the general formulas: 





= Vt[^ 


-ut- 


- iux - 


- TJUy] + U2t + iUxt + VUyt 




= Vx[(t> 


'-ut- 


- iux - 


- rjUy] + utx + iu2x + rjUxy 


(py 


= T^y[^ 


'-ut- 


- iux - 


- rjUy] + Uty + iUxy + r]U2y 




= T^xy [< 


i>-ut 


- iux 


T^liy] -\- Utxy ~\~ ^Uxxy ~1~ "H^xyy 


px 


= 'D2x[< 


i>-ut 


- iux 


V^y] ~l~ '^txx ~l~ C^xxx ~l~ V^xxy 


?y 


= T^2y[' 


i>-ut 


- CUx 


/yWy] ~l~ '^tyy ~t~ ^'^xyy ~1~ V^yyy 



(17) 



By extending the relations ()17p . substituting them into the condition ()16p and then equating 
with zero the coefficient functions of various monomials in derivatives of u, the following partial 
differential system with 11 equations is obtained: 

= C 

= Vu 

= Br]^ - D(f)2u 
= BCy-C<P2u 

= Arjy - TjAy - Au(t) + A^^ - ^Ax + +2D^y + 2Cri^ - At 

= Ari^ + 2Dr]y-riDy-^Dx-Du^-Dt 

= -A(t>2u + Bix-B(t>u + Briy-Bt-Bxi-Bu(t>-Byri 

= -T]t + Fr,,-Bct>x + Er,y-Et-Ex^-Eyrj-Eu(l) (18) 

+AVxy - Mxu + C^x + Dri2y " '^D(f>y^ 

= -(t-B(Py + F^, + ECy-Ft-FxC-Fyr]-F^(P 

^Cxy ~ ^4>yu + Ci2x + ^^2y ~ ^C'^xu 

= ct>t + G^^-Fct>x-E<l)y-Gt-GxS.-Gyr,-Gu(t> 

-Mxy - C'<^2x - D(l)2y 

The number of equations and of unknown functions which appear in the system (jlSp is relatively 
high. Two approaches are now possible: (i) to find the symmetries of a given evolutionary equa- 
tion, which means to choose concrete forms for A{x,y,t,u), B{x,y,t,u), G{x,y,t,u), D{x,y,t,u), 
E{x,y,t,u), F{x,y,t,u), G{x,y,t,u) and to use the system (fT8]l in order to find the coefficient 
functions ^{x, y, t), rj{x, y, t) and (/)(x, y, t, u) of the Lie operator; {ii) to solve the system (fT8|) taking 
as unknown variables A{x,y,t,u), B{x,y,t,u), G{x,y,t,u), D{x,y,t,u), E{x,y,t,u), F{x,y,t,u), 
G{x,y,t,u) and imposing a concrete form of the symmetry group. The first approach represents 
the direct symmetry problem and it is the usual one followed in the study of the Lie symmetries 
of a given dynamical system. The second approach, (ii), represents the inverse symmetry problem 
and it is more special, allowing us to determine all the equations which are equivalent from the 
point of view of the symmetry group they do admit. 
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2 Applications 



In the next considerations, we will solve the direct and inverse Lie symmetry problems for two 2D 
nonlinear models: the Ricci flow model and the convective-difl'usion equation. 

2.1 The Lie symmetry problems for 2D Ricci flow model 

One of the most fruitful models used in study of the black holes and in the attempt of obtaining a 
quantum theory of gravity is connected with the Ricci flow equations [17], 

We will investigate a 2D model for the Ricci flow equation, a nonlinear parabolic equation 
obtained when the components of the metric tensor g^p are deformed following the equation: 

d 

Q^9al3 = -Ra(3 (19) 

where Rai3 is the Ricci tensor for the n-dimensional Riemann space. The metric tensor of the space 
gai3 will be connected with the Riemann metric in the conformal gauge: 

ds^ = Qa^dx'^dx'^ = i exp{$(X, Y, t)}{dX^ + dY^) (20) 

The "potential" ^(X,Y,t) satisfles the equation: 

|.*=A* (21) 

It has been noticed [18] that the equation ()2ip is pretty similar with the Toda equation describing 
the integrable interaction of a collection of two dimensional flelds {^i,i = 1, 2} coupled by a Cartan 
matrix (Kij): 

^i^,,e*^(^'^) = A<I>,(X,y) (22) 
j 

Introducing the held u{x,y,t) given by 

n(x,y,t)=e* (23) 

the equation ()2ip takes the form: 

ut = {lnu)xy (24) 

An equivalent form for the previous equation, which will be used in the next considerations of the 
paper, is: 

"ixy ^x^y 
1? 

The previous equation could be derived from the general one ()12p by choosing the following par- 
ticular coefficient functions: 

A{x, y, t, u) = -,B{x, y, t, u) = — \, 

C{x,y,t,u) = D{x,y,t,u) = E{x,y,t,u) = F{x,y,t,u) = G{x,y,t,u) = (26) 
2.1.1 Lie symmetries for 2D Ricci Row model 

To finding the Lie symmetry operators for the Ricci flow model ()25p we have to solved the general 
determining system (jlSp in the conditions (j26p . The solution is represented by the coefficient 
functions S,{x,y,t), rj{x,y,t), (f){x,y,t,u) which determine the Lie symmetry operator ()13p . It has 
the form: 



= ^ - ^ (25) 
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As U contains coefficients in the form of two arbitrary functions rj}, we deal with an infinite 
number of symmetry operators. The action of U can be split in various "sectors", depending on 
the concrete form we might choose for these functions. 

Let us consider the linear sector of the Lie symmetries in which the forms of the coefficient 
functions of the symmetry generator ()13p are: 

(p = I, C = "f^^ + ci, r] = vy + C2, 4> = — (m + v)u (28) 

with m, V, k, ci, C2, C3 arbitrary constants. 

The general Lie symmetry operator ()13p becomes: 



d d d d 

U{x,y,t,u) = + + ci)^ + {vy + ^2)^ - (m + v)u— (29) 

Consequently, the nonlinear Ricci flow equation ()25p admits the 4— dimensional Lie algebra spanned 
by the independent operators shown below: 

d d d d d d 

Vi=x^-n^,V2 = ^, V3 = y^- u^, V, = ^ (30) 
ox ou ox oy ou oy 

The forms of the operators Vi, i = 1,4 suggest their significations: V2,Va generate the symmetry 
of space translations, Vi, C/3 are associated with the scaling transformations. 

When the Lie algebra of these operators is computed, the only non-vanishing relations are: 

[V2,Vi] = V2AyA,V^] = Vi (31) 
2.1.2 Optimal system of subalgebras for 2D Ricci fiow^ model 

It is well known that reduction of he independent variables by one is possible using any linear 
combination of the generators of symmetry (j30p Vi,i = 1,4. We will construct a set of minimal 
combinations known as optimal system [5j. To construct the optimal system we need the commu- 
tators of the admitted symmetries given in the Table 1. 

[VuVj\ Vi V2 ^3 ^4 

Vi -V'2 

^2 V2 Tablel: Lie brackets of the admitted symmetry algebra 

-Vi 
Vi ^4 
An optimal system of a Lie algebra is a set of Z— dimensional subalgebras such that every 
/—dimensional subalgebra is equivalent to a unique element of the set under some element of the 
adjoint representation. The adjoint representation of a Lie algebra {Vi,i = 1,...,4} is constructed 
using the formula [5J: 

n 2 

Ad{eMeV^))V, = ^{o^dVrVj = V, - e[Vu V,-] + -[V, [V, V,]] - ... (32) 

n 

Let us consider the linear combination of the symmetry generators: 

^ = 01^1-^02^2 + 03^^3 + 04^4 (33) 

Our task is to simplify as many of the coefficients Oj as possible through judicious applications of 
adjoint maps to ^.Suppose first that oi / in (f33|l . One may re-scale oi such that ai = 1. We 
start with the combination: 

= Vi + a2V2 + 031^3 + 041^4 (34) 
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If we act on we can make the coefficient of V2 vanish: 

y(2) = 14 + asVs + a^Vi (35) 
Next, we act on F^^) by Ad{eyL\){^Vij) to cancel the coefficient of V4, leading to the operator: 

y(3) = v^ + 03^3 (36) 



Using the adjoint representation (|32|1 . no further simplification is possible. Consequently, the 
1— dimensional subalgebra spanned by V with oi 7^ is equivalent to the one spanned by Vi + /3V3, 
/3 G i?. 

The remaining 1— dimensional subalgebras are spanned by operators with oi = which have 
the expressions: 

= a2V2 + aaVg + 04^4 (37) 

Let us assume that 02 7^ and scale to make 02 = 1. Now we act on by Ad{Qiq){^Vi)) so that 
it is equivalent with the operator: 

= V2 + a^V^ (38) 

No further simplification is possible. Consequently, the 1— dimensional subalgebra spanned by V 
with a2 7^ is equivalent to the one spanned by V2 + aVa, a £ R. 

If we consider the case ai = 02 = 03 = 0, 03 / 0, 03 = 1, the following generator is obtained: 

= V3 + a^Vi (39) 

Acting on V^^^ by ^d(exp(a4V4)), we obtain the operator V3 which represents the next subalgebra 
of the optimal system. 

Finally, let us consider the last case oi = 02 = 03 = 0, 04 7^ 0,04 = 1 in (f33|) . Results the last 
subalgebra V4. 

In conclusion, the optimal system of 1— dimensional subalgebras has the form: 

{y2 + a^3,^i+/3^3,^3,^4} (40) 

2.1.3 Invariant solutions for 2D Ricci flow^ 

Let us pass now to the problem of the invariant quantities. We shall analyze the invariants 
associated with the optimal system of symmetry operators (|40p. 



(41) 



• The operator V2 + aVs from ()40p has the characteristic equations: 

dt dx dy du 



1 ay —au 
By integrating these equations result 3 invariants with expressions: 

h = t, l2 = ye-"^ /3 = yu (42) 

By introducing the similarity variable z = I2 = ye~"^, designating the invariant /3 = h{t,z) 
as a function of the other ones, the following solution is obtained: 

uit,x,y) = ^ (43) 

Setting the derivatives of ([33]) into the Ricci equation ([25]) . results the similarity reduced 
equation for h{t, z) with the form: 

hth^ — az^hh2z — az'^hj, + azhhz = (44) 
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The solution of the previous equation is: 



with a,ri,r2,r^ arbitrary constants and z the similarity variable. 
Consequently, the invariant solution corresponding to operator V2 + aVs has the final form: 

u{t, x,y) = -- [r,t + + tanh ^ ^- J J (46) 

• The operator Vi + /SVJj from (j4Up has the characteristic equations: 

dt dx dy du 



X -(l + /3)n 



(47) 



In this second case, following the same procedure, we obtain also 3 independent invariants 
with expressions: 

h = t, l2 = yx-^ h = y^'^f'^/^u (48) 

By introducing the similarity variable z = I2 = yx~^, designating the invariant = h{t,z) 
as a function of the other ones, the following solution is obtained: 

u{t,x,y) = h{t,z)y-'^^+'^^/f^ (49) 

Setting the derivatives of (|I9|) into the Ricci equation ([25]) . results the following (1 + 1) reduced 
equation for h{t,z): 

hth^z^-^/^-^'^ + I3hh2z + f3z-^hh, -I3hl = (50) 
The solution of the previous equation is: 

with f3,pi,p2,P3,P4: arbitrary constants and z the similarity variable. 

Consequently, by using the invariant solution corresponding to the operator Vi + PV^ has 
the final form: 



Because ()25p is symmetric in x and y, there is also a second similarity solution of the form: 

u{x,y) = ^^M^ ygsix) (53) 

y 

which is generated by the symmetry operator V3 from (|3U|) . 

Again, by the reason of symmetry in x and y of the analyzed model (j25p . the last similarity 
solution associated to the symmetry operator V4 from ([30]) . is generated as below: 



u{x) = Qiix), 'igA{x) (54) 
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2.1.4 The class of equations with Ricci type symmetries 



Now, our aim is to find the class of equations with generic form (|12p which admit Ricci type symme- 
tries ()28p . Consequently, we have to solve the system (jlSp taking as unknown functions A(x, y, t, u), 
B{x, y, t, u), C{x, y, t, n), D{x, y, t, u), E{x, y, t, u), F{x, y, t, u), G{x, y, t, u) and imposing a concrete 
form of the symmetry group. 

Let us consider the linear sector for the Lie symmetries where coefficient functions of symmetry 
operators have the general expressions: 



(p = 1, C = mx + ci, r] = vy + C2, (f) = ku + C3 



(55) 



with m, V, k, ci, C2, C3 arbitrary constants. 

Remark 1: For the concrete values k = — (m + v),C3 = and nonvanishing values for all 
the other constants m,v,k,ci,C2, the linearization ([55]) is reduced to (j28|) considered for the Ricci 
flow model. The results will be particularized for this case when the Ricci symmetries have been 
imposed. 

With the choice (|55p . the general differential system (|18p has for the "unknown" functions 
A{x,y,t,u), B{x,y,t,u), C{x,y,t,u), D{x,y,t,u), E{x,y,t,u), F{x,y,t,u), G{x,y,t,u) solutions 
which can be expressed in terms of 14 arbitrary functions {J'j, ^j^j = They are the 

general forms: 



A{x,y,t,u) 

B{x,y,t,u) 
C{x,y,t,u) 

D{x,y,t,u) 
E{x,y,t,u) 
F{x,y,t,u) 
G{x,y,t,u) 



T[ +j^f(x",y",z")??— 
T'2 [X', Y\ z') + -^'(X" , y" , z")7? 

J-^ {X', Y', Z') e + -^'(X" , Y^' , ^" w 
T'^ {X' X ,Z')i^ ^ Tl{X\Y\Z'')rf^ 
T'^ [X', Y\ Z') + J^'(X" , y" , Z" )r/ V 
J-^ {X', Y', Z') + J^'(X", y", Z")r/v 
-F7 {X',Y',Z')C^ + J-7'(X",y",Z")r/^ 



m-\-v — k 



(56) 



where 



X' = 



X" = 



""^ y, ^ "rat - Ing ^, 



m 



(f)^{ m ) 
k 



m V k 

Remark 2: In the case A; = 0, the arguments Z', Z" are changing and they become: 

^, ^ mu - C3 In g ^„ ^ vu - C3 In r? 
m ' V 

It is useful to choose particular solutions. Let us consider the cases: 
Casel : C3 = and (m + v) is a multiple of the constant k, that is: 

m + v = k-n, (V)n G Z\{0} 



(57) 



(58) 



(59) 



Following the general solution ([56]) - ([57|) . the concrete forms for A{x, y, t, u), B{x, y, t, u), C{x, y, t, u), 
D{x,y,t,u), E{x,y,t,u), F{x,y,t,u), G{x,y,t,u) may be in this case: 



A{x,y,t,u) 
B{x,y,t,u) 



2 ^ 



2 ' 



Ik) 



+ 



(n-l) 



(60) 
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C{x,y,t,u) = D{x,y,t,u) = E{x,y,t,u) = F{x,y,t,u) = G{x,y,t,u) = 



where the parameter k has the form (j59p . 

The 2D general equation (|12p takes the particular form: 



(61) 



Remark 3: For n = — 1, the previous equation generates the 2D Ricci flow model ()25p . 
Case2 : k = 0, = m + v and nonvanishing values for all the other constants m, v, ci, C2. 
The solution (1561) could take in this case the form: 



A{x,y,t,u) = 


2 ^ 




2 ' 


B{x,y,t,u) = 


2 ^ 




2 ' 



{-) 



(^) = — = 
du 



(62) 



C(x, y, t, n) = D{x, y, t, u) = E{x, y, t, u) = F{x, y, t, u) = G{x, y, t, u) = 
In this second case, the following 2D evolution equation of type (|12p is generated: 

Ut — 6 Uxy ~t~ 6 Ux^y — (^ ^x)y 



(63) 



Remark 4: The equations ([6T]1 and ([63]) generated by the previous two cases correspond to 
the 2D nonlinear heat equation which has the general expression 



Ut = {g{u)ux)y 



(64) 



It was proven [21] that the choices g{u) = and g{u) = are the only two possible cases for 
which the Lie symmetries exist and more the Lie operators have linear forms. 



2.2 The Lie symmetry problems for 2D convective-diffusion equation 

The second nonlinear application is represented by the 2D convective-diffusion equation [?]. It is 
a parabolic partial differential equation, which describes physical phenomena where particles or 
energy (or other physical quantities) are transferred inside a physical system due to two processes: 
diffusion and convection. In the simpler case when the diffusion coefficient is variable, the convection 
velocity is constant and there are no sources or sinks, the equation takes the form: 



Ut = UU2x + UU2y — VUx (65) 

with diffusion coefficient u and convective velocity v = const, belongs to the Ox direction. 

It is easy to remark that (j65p results from the general class of equations (jl2p by choosing the 
particular functions: 

C(x, y, t, u) = D{x, y, t, u) = u, F{x, y, t, u) = -v 

A{x,y,t,u) = B{x,y,t,u)=E{x,y,t,u) = G{x,y,t,u) = (66) 
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2.2.1 Lie symmetries for 2D convective-diffusion equation 



In the conditions (|66p the general determining system (jl8|) for symmetries becomes: 

<p2u = 

2u^,-(t) = 

2ur]y-<i) = (67) 
- + U7?2x- + uri2y - 2u(t)y^ = 

It has the solution: 

^ = —{x - vt) + C2y + C3, r] = —y - C2{x - vt) + C4, 4> = ciu (68) 

In this case, the Lie symmetry generator takes the expression: 

rr/ N 9 (C\, , \ d ( c\ , . \ d d , . 

U{x,y,t,u) = -qI + y-i^{x - vt) + C2y + c^j + y—y-C2{x-vt) + C4j ^ + ^1^^ (^9) 



Consequently, the nonlinear convective-diffusion equation ()65p admits the 4— dimensional Lie alge- 
bra spanned by the operators shown below: 

fx — vt\dfy\d d 

= + (I) % + <™> 

d , . d d d 

V2 = y—-{x-vt) — ,V:i = —,Vi = — 
ox oy ox oy 

When the Lie algebra of these operators is computed, the only non-vanishing relations are: 

[V^, Vi] = \vz, [Va, Vi] = ^4, [V2, Vz] = ^4, [^4, V2] = (71) 
2.2.2 Optimal system for convective-diffusion equation 

For this model the commutators of the symmetry operators ()70p are given below in the Table 2: 



Table 2: Lie brackets of the admitted symmetry 



algebra 

Let us consider the linear combination of the symmetry generators: 

V = hiVi + 62^2 + ^3^3 + ^4^4 (72) 

Our task is to simplify as many of the coefficients 6j as possible through judicious applications of 
adjoint maps to ^/.Suppose first that 61 7^ in ([72]) . One may re-scale hi such that hi = 1. We start 
with the combination: 

= Vi + 62^2 + 63^3 + ^4^4 (73) 





Vi 


V2 


^3 


V4 


Vi 








-^3/2 


V4 


V2 










-^3 




1^3/2 














^3 
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If we act on we can make the coefficient of V3 vanish and we obtain the 

operator: 

y{2) =Vi + b2V2 + b'^Vi, b'^ = bi + 26263 (74) 

Using the adjoint representation ()32p for our model, no further simplification is possible. Conse- 
quently, the 1— dimensional subalgebra spanned by V with 61 7^ is equivalent to the one spanned 
by Vi + aV2 + V a, /3 G i?. 

The remaining 1— dimensional subalgebras are spanned by operators with 61 = which have 
the expressions: 

= b2V2 + 631/3 + 6414 (75) 
Let us assume that 62 7^ and scale to make 62 = 1. Now we act on 

by Ad{exp{b3V4)) so that 

it is equivalent with the operator: 

=V2 + 6414 (76) 

Here, further simplification is possible. If we act on by ^d(exp(-64y3))- In this case, we 
obtain the operator V2 which is the following 1— dimensional subalgebra spanned V with 61 = 
and 62 = 1. 

If we consider the case 61 = 62 = 0, 63 / 0, 63 = 1, the following generator is obtained: 

y(3) = ^3 + 64-14 (77) 

If we act on V^^^ by Ad{exp{eV2)) , where e is the solution of the equation 

je^+e-b4 = (78) 

we can vanish the coefficient of V4 and we obtain the operator: 

V^^^=b'3V3, 6'3 = 1 + 646- y (79) 

where e verified (|78p . 

Consequently, by setting 61 = 62 = 0, 63 = 1 in ([72]) is generated V3 which is the last subalgebra 
of the optimal system. 

In conclusion, the optimal system of 1— dimensional subalgebras for 2D convective-diffusion 
equation is: 

{V2,V3,Vi + aV2 + pV4,ya,(3€R} (80) 



2.2.3 Invariant solutions for the convective-diffusion equation 

Through the reduced similarity method, each operator {V^, i = 1, 4} can generate invariant solutions 
of the model. Let us illustrate for our case what are the concrete forms of the similarity solutions 
generated not by this base of operators, but by the set of the optimal ID subalgebras (|8Up . 



(81) 



Taking into account ()70p . the symmetry operator V2 has the characteristic equations: 

dt dx dy du 
y vt — X 



In this second case, following the same procedure, we obtain also 3 independent invariants 
with expressions: 

h=t, l2 = vtx-^-^, l3 = U (82) 
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With the notation I2 = z and = u = g{t, z), the reduced equation for g{t, z) wih take the 
form: 

gt + {2z - vh^)gg2z + 2gg, + vhg, = (83) 

It admits the solution: 

2z - vH"^ + 2qi 



where qi,q2,v are arbitrary constants. 



4i + 2q2 



(84) 



Thereby, the second similarity solution corresponding to operator V2 has the final form: 

2vtx — — y"^ — v^t^ + 2qi 

u(t,x,y) = : 

^ 4i + 2q2 

• The operator V3 from ()70p yields the characteristic equations: 

dt dx dy du 

0" ~ T ~ T ~ T 

Therefore also 3 invariants are generated: 

h = t, l2=y, h = u 



(85) 



(86) 



(87) 



Once again, expressing the last invariant I3 as a function of the others ones, we obtain the third 
similarity solution: 



u{t,y) 



+ l^y + 94 
92 - qit 



(88) 



with 91,92 arbitrary constants. 



Again on the basis of ([70]) . the last operator from ([80]) . Vi + aV2 + f3Vi, has the characteristic 
equations: 

dy 



dt 



dx 



du 



ay + ^ \ + a{vt - x) ^ ^ u 
By integrating these equations one obtains 3 invariants with expressions: 

I + a{vt - x) + /3 _ 



h = t, h 



u 



\ + a{vt — x) + 13 



(90) 



By introducing the similarity variable z = I2, designating the invariant = h{t, 
function of the other ones, the following solution is obtained: 



u{t,x,y) = h{t,z) 



y + a{vt - x) + /3 



n 2 



(91) 



Setting the derivatives of (|9T|) into the convective-diffusion equation (|65|) . we obtain the 
following (1 + 1) reduced equation for h{t, z) : 



ht-2{a^ + ^] z^hhz - 4 ( + t 1 zhh^ -2{a^ + ^]h^ =0 



1 



1 



The solution of the previous equation is: 

h{t,z) = 



2(a2 + i)t-7 



(92) 



(93) 
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with a, 7 arbitrary constants. 

Using ([9T]l . the invariant sohition generated by the operator Vi + aV2 + /3Vi is pointed out: 

-I 2 



u{t,x,y) 



1 



y 



+ a{vt - x) + /? 



2(a2 + i)t-7 
where a,/3,7,t; arbitrary constants. 
2.2.4 Inverse symmetry problem for 2D convective-diffusion equation 



(94) 



Our aim is now to find the class of equations with generic form (jl2p which admits the same symme- 
tries with those corresponding to 2D nonlinear convective-diffusion equation ()65p . Consequently, we 
have to impose that the coefficient functions (j68p which determine the base of symmetry operators 
()70p verify the general determining system (jlSp . 

The solutions of differential system (jlSp describe the coefficient functions of the general evolu- 
tionary equation (fl^ as follows: 



A 


= B = 


0, c = 


L> 


E{u) 


= a/u 


C4 COS ^ 


C2 








F{u) 


= ^/u 


C4 sin ^ 


C2 




Cl 


G{u) 


= CqU 







C2 
Cl 
C2 
Cl 



(95) 



where Cj, j = 1,6 and u are arbitrary constants. 

In particular, for C3 = 1, C4 = C5 = cg = and arbitrary ci and C2, the solution ([95]) generates 
the 2D nonlinear convective-diffusion equation ()65p discussed above. 



3 Conclusions 

This lecture intended to present some key aspects on how a dynamical systems whose evolution 
is described by a nonlinear differential equation can be studied using the symmetry method. The 
main steps which have to be done in order to find a set of exact solutions are: (i) determination of 
the general form for the symmetry operator; (ii) determination of the optimal set of independent 
operators which can generates the minimal subalgebras; {in) based on the optimal set of indepen- 
dent operators and using the similarity reduction procedure, a complete set of invariant solutions 
can be generated; {iv) last but not least, a special method can be applied in order to find the largest 
class of nonlinear differential equations which belong to the same class as a given equation in the 
sense of the symmetries they observe. This algorithm was applied for two important examples of 
nonlinear 2D partial derivative equations, the Ricci flow and the convective-diffusion equation. For 
the first example the optimal system of subalgebras contains the same number of generators, four, 
as the whole symmetry algebra. The optimal system of symmetry subalgebra for the convective- 
diffusion equation has the dimension three, despite the existence of four independent symmetry 
operators. In both cases, the whole set of invariant solutions had been obtained. 
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